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ABSTRACT 

GRB 080503 is a short gamma-ray burst (GRB) detected by Swift and has been classified as a compact-star- 
merger-origin GRB. The soft extended emission and the simultaneous late re-brightening in both the X-ray and 
optical afterglow lightcurves raise interesting questions regarding its physical origin. We show that the broad¬ 
band data of GRB 080503 can be well explained within the framework of the double neutron star merger model, 
provided that the merger remnant is a rapidly-rotating massive neutron star with an extremely high magnetic 
field (i.e. a millisecond magnetar). We show that the late optical re-brightening is consistent with the emission 
from a magnetar-powered “merger-nova”. This adds one more case to the growing sample of merger-novae 
associated with short GRBs. The soft extended emission and the late X-ray excess emission are well connected 
through a magnetar dipole spin-down luminosity evolution function, suggesting that direct magnetic dissipation 
is the mechanism to produce these X-rays. The X-ray emission initially leaks from a hole in the merger ejecta 
pierced by the short GRB jet. The hole subsequently closes after the magnetar spins down and the magnetic 
pressure drops below ram pressure. The X-ray photons are then trapped behind the mergernova ejecta until the 
ejecta becomes optically thin at a later time. This explains the essentially simultaneous re-brightening in both 
the optical and X-ray lightcurves. Within this model, future gravitational wave sources could be associated 
with a bright X-ray counterpart along with the mergernova, even if the short GRB jet beams away from Earth. 

Subject headings: gamma rays: burst - hydrodynamics - radiation mechanisms: non-thermal - stars: neutron 


1. INTRODUCTION 

The frequency range of the next generation gravitational- 
wave (GW) detectors, suc h as Advanced LIGO (lAbbott et alJ 
1200^, Advanced VIRGO (lAcernese et al.l2008h and KAGRA 
(iKuroda & LCGT Collaborationl20I(ill interferometers, is de¬ 
signed to uncover the final inspiral and merger of compact 
object binaries (NS-NS, NS-BH, BH-BH systems). Due to 
the faint nature of GW signals, an associated electromagnetic 
(EM) emission signal coinciding with a GW signal in both 
trigger time and direction could play a crucial role for con¬ 
firming the astrophysical origin of the GW signals and study¬ 
ing the astrophysical origin of the GW sources (e.g. host 
galaxy, distance, etc). 

Short-duration 7-ray bursts (SGRBs) have long been pro- 
pose d to originate fr o m mergers of cornpact object bina¬ 
ries (iPac zvnsri 11 9861 : lEichler et alJ Il989t iPaczvnskil 119911 : 
iNaravan et al.lll993f If so, SGRBs may provide the bright¬ 
est EM counterpart associated with events detected by those 
upcoming interferometers. However, observations of SGRBs 
suggest that at least some of them are collimated into a smal l 
opening angle (iBurrows et alJl2()()^lDe Pasquale et al.ll20ir)l) . 
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SO that most GW signals would not b e detected together with 
SGRBs (e.g. lMetzger & Bergeill201^ . Lately, additional EM 
signatures of the compact binary mergers (espe cially for NS- 
NS system) becomes a topic of growing interest (lBerged20l4 
for a review). 

Numerical simulations show that a mildly isotropic, sub- 
relativistic outflow could be ejected during the merger of 
binary neutron stars, including the tidal tail matter during 
the merger and the matter from the acc r etion disk (e.g. 
iRezzolla et al.l 12 01 It [Rosswog et akl 1201 3t iBauswein et alJ 
2013t[Hotokezaka et alJ2013l) . The typical mass and speed of 
the ejecta ar e in the range of 10~^- 10~^ Mq and 0.1 -0.3 c, 
respectively (iHotokezaka et al.ll2013h . Recently, several inter¬ 
esting EM signatures from the ejecta have been well studied, 
whose brightness are essentially determined by the properties 
of the left over remnant from the merger. 

Usually, the merger product is assumed to be either a 
black hole or a temporal hyper-massive neutron star which 
survi ves 10-100 ms befor e collapsing into t h e black hole 
(e.g. iRosswog et akl 120031: lAlov et all l2005l IShibata et alJ 
l20()5t IRezzolla et al.l 1201 It IRosswog et al.1 l2013h . In this 
case, an optical/infrared transient is expected to be pow- 
ered by radioactive decay from r-process radioactive rnaterial 
jLi_&Pacz^nsy 19^ lKulkarnill2005t iMetzger et alJ 120101 : 
[Barnes & KasenlT2013h . henceforth we call it a r-process- 
powered merger-nova^. Besides this thermal emission, a 
long-lasting radio emission is also expected from the in¬ 
teraction between the ejecta and the ambie nt medium, al¬ 
though it is normally too weak to be detected (iNakar & PiranI 
I2OI It IMetzger & Bergei]l20T^lPiran et al.ll2013h . Such tran- 
sients are more isotropic than SGRBs. Depending on the 
direction of our line of sight, these transients could be de¬ 
tected alone or to be accompanied by SGRBs, provided 
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that t heir luminosities are large enough ([Metzger & Bergeil 
2Q12h. After several years of search (iBloom et alJ |200 ct 
P erlev et al.ll2009l : lKocevski et al.l2010l) . an r-process-powerd 
merger-nova was finally claimed to be detected in the in¬ 
frared band with the Hubble S pace Telescope (HS T) for GRB 
1 ^OhO ^B dTanvir et al.llpi Berger et alTmT^ . More re¬ 
cently, I Yang et alJ (1201 5h re-examined the late afterglow data 
of GRB 060614 observed with HST, and found a significant 
F814W-band excess at f ^ 13.6 days after the burst. They 
claimed that it is very likely another candidate of r-process- 
powered merger-nova. For both cases, the mergernova inter¬ 
pretation was based on one single data point. 

Alternatively, it has long been proposed that the post¬ 
merger product could be a stable or super-massive mil¬ 
lisecond magneter, if the equation of state of nuclear mat¬ 
ter is stiff enough and t he total mass of the two neu- 
tron stars is smal l enoug h (iDai et alJ 12005 iFan & Xu|l200^ 
iGao & F^ 120061: IZhanal^Ol : iGiacomazzo & Pernall2013h . 


Evidence of a magnetar following some SGRBs has been 
collected in the Swift da t a, including the exte nded emission 


Norris & Bonnelll 120061 iMetzser et al. 

I2OO8I), X-ray flares 

(Barthelmv et alj|2005l ICamnana et al.l 

20061) and more im- 

portantly, “internal plateaus" with rapid decay at the end of 


the plateaus (iRowlinson et alJ 1201 Ol 12013L iLii et alJ l2015h . 


Nonetheless, available observations (e.g., the lower limit of 
the maximum mass of Galactic NSs and the total mass dis¬ 
tribution of Galactic NS-NS binaries) and numerical simu¬ 
lations al low the existence of the post-merg er massive NS 
remnant dZhangl 120131 : iMetzger & Ph^ 120141 and reference 
therein). Compared with the black hole merger remnants, the 
main consequences of magnetar merger remnants include the 
following: 

• The spin-down of the NS remnant supplies an addi¬ 
tional energy source to the system; 

• The strong neutrino-driven wind from the NS provides 
additional mass outflow to the system. 

In this case, the detectable EM signatures from the system be¬ 
come much richer. Besides the putative short GRB signature, 
the following EM signals may be expected. Eirst, the mag¬ 
netar would eject a near-isotropic Poynting-flux-dominated 
outflow, the dissi pation of whi ch could power a bright early 
X-ray afterglow (IZhang|l2013h . Second, the thermal emis¬ 
sion from the ejecta could be significantly enhanced due 
to additional heating from magnetar wind (lYu et al.l 120131 : 
iMetzger & Pirol2014l) . This power could exceed the r-process 
power, so that we may call the corresponding transient as 
a “magnetar-powered merger-nova”. Einally, the magnetar- 
power would energize and accelerate the ejecta to a mildly 
or even moderately relativistic speed, and the interaction be¬ 
tween the ejecta and the ambient medium could produce a 
strong external shock that gives rise to bright broad-band 
emissi on (i.e. the double neutron star (DNS) merger afterglow 
model. iGao et al.l2013a^ . Some recently discovered transients 
could be interpreted within such a scenario, lending support 
to a post-merger magnetar remnant. Eor instance, the Palo- 
mar Transient Eactory (PTE) team recently reported the dis¬ 
covery of a rapidly fading optical transient source, PTEl lagg. 
Lacking a high-energy counterpart, it has been proposed to 
be a good candi date for the DNS merger afterglow emission 
(IWu et alj|2014l) . Moreover, considering its broad-band data, 
GRB 130603B and its claimed “kilonova” can be interpreted 
within the framwork of a magnetar-powered DNS merger 


remnant given that the mag netar underwent s i gnificant energy 
lost th rough GW radiation (lEan et alJl2013bl: Metzger & Pird 
12014 . 

Similar to GRB 130603B, GRB 080503 is a SGRB with 


bright extended emission. Based on its negligible spectral 
lag of prompt emission and extremely faint afterglow, GRB 
080503 has been classified as a compact-star-merger-origin 
GRB^° (iPerlev et al.l l2009ll . The most peculiar feature in 
GRB 080503 is that after the prompt emission (began with a 
short spike and followed by extended emission) and the early 
steep decay afterglow phase, it didn’t immediately enter into 
the regular afterglow phase. Being signal-less for about one 
day, it presented a surprising re-brightening in both the op¬ 
tical and X-ray bands. In the optical, it remained bright for 
nearly five days. Within the post-merger remnant is a black 
hole, the scenario has been investigated for GRB 080503. A 
“r-process-powered merger-nova" model can marginally ex¬ 
plain t he optical data, but the X-ray data coul d not be inter¬ 
preted (iPerlev et alJl2009l : iHascoet et al.ll2012h . In this work, 
we make a comprehensive analysis on the multi-band obser¬ 
vations of GRB 080503, and suggest that the magnetar merger 
remnant scenario can well interpret the entire data set, making 
a solid case to connect the late optical excess of GRB 080503 
with a magnetar-powered merger-nova. We note that the idea 
that GRB 080503 is a good candidate for a m agnetar-powered 
transie nt has been qualitatively proposed by IMetzger & Pirol 
(l20T4 . 


2. OBERVATIONAL FEATURES OF GRB 080503 

GRB 080503 was detected by Burst Alert Telescope (BAT) 
on abroad Swift satellite at 12:26:13 on 2008 May 3 (see ob¬ 
servational details in IPerlev et al.l (l2009l) i. Its prompt emis¬ 
sion (in the 15-150 keV bandpass) contains a short bright ini¬ 
tial spike with a duration of 0.32 ± 0.07 s, followed by a soft 
extended emission lasting for 232 s. The peak flux of the 
initial spike (measured in a 484 ms time window) and the flu- 
ence of the extended emission (measured from 5 s to 140 s 
after the BAT trigger) are (1.2 ± 0.2) x 10“^ erg cm“^ s“' and 
(1.86 ±0.14) X 10“® erg cm respectively. Although the 
fluence ratio between the extend ed emission and the spike is 
as large as 30 in 15- 150 keV, IPerlev el'iD (l2009l) further 
analyzed its features of hardness ratio and spectral lag in de¬ 
tail, and found that this burst is still more reminiscent of a 
compact-star-merger-origin GRB. 

After the extended emission phase, the X-ray light curve de¬ 
cays rapidly (a = 2-A, where Tv ex; f““) until below the XRT 
detection threshold, and kept undetectable for about 1 day (as 
shown in Eigure|2|l. Then the X-ray flux rebrightened to the 
level of 10“^ /iJy around 10^ s after the BAT trigger. 20 days 
later, the Chandra X-Ray Observatory ACIS-S was again em¬ 
ployed to conduct imaging on the relevant position, but no 
source was detected. 

In the optical band, many facilities were employed to search 
for afterglow signals on the first night after the trigger, such as 
Swift UV-Optical Telescope (UVOT), Keck-I telescope and 
Gemini-North telescope, only a single Gemini g band de¬ 
tection was obtained at 0.05 day. However, on the second 
night after the trigger, the afterglow surprisingly rise above 
the detection threshold to the level of 10“^ /iJy and kept 
bright for nearly five days. Later on, the localization region 

The physical category of a GRB may not always be straightforwardly 
inferred based on the duration information, and multi-band observational cri¬ 
teria are needed (Zhang et al. 2009). 
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was observed with Hubble Space Telescope in two epochs on 
2008 May 12 and July 29. Although only upper limits were 
achieved, the results infer a rapid decay feature for the late 
optical excess component. 

During the observation, many attempts to measure the red- 
shift of GRB 080503 were operated, even with Hubble Space 
Telescope. Unfortunately, only an upper limit, z < 4, was 
achieved. 


3. MODEL DESCRIPTION 
3.1. General picture 

If the equation of state of nuclear matter is stiff enough, 
the central product for a binary neutron star merger could 
be a stable or a supra-massive NS rather than a black hole. 
This newborn massive NS would be rotating with a rota¬ 
tion period in the order of milliseconds (close to the cen¬ 
trifugal break-up limit), and may also contain a strong mag¬ 
netic field B > 10*‘*G similar to “magnetars" (IZhandl2013l: 
iMetzger & Pirol20l4 and reference therein). The millisecond 
magnetar is surround by a sub-rela tivistic (vej ^ 0.1-0.3c ) 
ejecta with mass ^ (10“"^- 1O“^)M0 (iHotokezaka et al.l2013l) . 
Considering a variety of the origins for the ejecta materials, 
a sphe rical symmetry could be reasonably assumed for the 
ejecta (iMetzger & Piroll2014l) . 

Shortly after the formation, the magnetar would be 
surrounded by a c e ntrifugally supporte d accretion disc 
dMetzger et al.l 120081: iDessiut et al.l 120091: iLee et alJ 120091: 
iFernandez & Metzgeiil2013b, launch i ng a s hort-li ved (< s) 
collimated iet (IZhang & E)Sl l2008l l2009l l2010h . The jet 
could easily punch through the ejecta shell and then power 
the prompt short sp ike emission and the broad band GRB 
afterglow emission (IMetzger et al.l l200^ . After the whole 
jet passing through, it is possible that the gap remains open 
as the Poynting-flux-dominated magnetar wind continuously 
penetrates through the hole. Due to the dynamical motion 
of the ejecta, the ejecta materials tend to quench the out¬ 
flow by closing the gap. During the early spindown phase 
when the Poynting-flux luminosity is essentially a constant, 
the ram pressure of the ejecta may be balanced by the mag¬ 
netic pressure of the outflow. After the characteristic spin- 
down time scale fsd, the magnetic pressure drops quickly, so 
that the gap is closed in a time scale of fciose- The total dura¬ 
tion when the magnetar wind leaks from the ejecta and make 
bright X-ray emission due to internal dissipation (e.g. through 
an internal-collision-induce d magnetic reconnec tion and tur¬ 
bulence (ICMART) process, IZhang & Yanll201 Tl) is the dura¬ 
tion of the extended emission, i.e. fee = tsd + Mose- According 
to a more detailed estimate (see Appendix for details), this 
duration can be consistent with the observed duration of ex¬ 
tended emission given reasonable parameters. After fee, the 
magnetar wind is stifled behind the ejecta, so that soft ex¬ 
tended emission stops and the high-latitute emission in X-ray 
band shows up, which gives a rapidly dropping X-r ay tail 
dKumar & Panaitescull2000l: IZhang et al.ll2007l l2009l) . It is 
similar to the situation when the central engine is shut down. 

Trapped by the ejecta materials, the magnetar continuously 
spin down, and when the magnetar wind encounters the ejecta, 
a signihcant fraction of the wind energy (parameterized as 0 
could be deposited into the ejecta, either via direct energy in¬ 
jection by a Poynting flux (iBucciantini et al.ll2012h . or due to 
heating from the bottom by the photons generated in a dissi¬ 
pating magnetar wind via for ced reconnect ion (if R < Rdis) or 
self-dissipation (if R > Rjis) (IZhang||2013l) . Such continuous 


energy injection n ot only heats the ejecta material to p ower 
the merger-nova, (lYu et al.ll2013l: IMetzger & Piroll2014l) . but 
also accelerates the ejecta to a mildly or moderately relativis¬ 
tic speed, giving rise to strong afterglow emission by driving 
a stron g forward shock into the ambient medium (iGao et al.l 
l2013ah . Note that the remaining fraction of the wind energy 
(1-0 would be stored in the trapped dissipation photons and 
eventually diffuse out with a deducing factor e~'^, where r is 
the optical depth of the ejecta. 

In summary, there are four emission sites and several emis¬ 
sion components involved in this model (as shown in Figure 
[TJ: i) the jet component that powers the short spike in prompt 
emission and the GRB afterglow emission; ii) the early mag¬ 
netar wind component that powers the soft extended emission 
and the high latitude tail emission; iii) the magnetar-powered 
merger-nova emission acomponent, and the corresponding 
DNS merger afterglow emission; iv) the late magnetar-wind- 
powered X-ray component when the ejecta becomes transpar¬ 
ent. In the following, we describe the details for calculating 
these main emission components. 

3.2. Magnetar wind dissipation 

Considering a millisecond magnetar with an initial spin pe¬ 
riod Pi and a dipolar magnetic held of strength B, its total 
rotational energy reads 0ot = (l/2)ml ~ 2 X 10 ^2/45^^2^ erg 
(with /45 ^1.5 for a massive neutron star). The spin-down 
luminosity of the magnetar as a function of time could be ex¬ 
pressed as 


/ t \ ^ 

Tsd — Tsd.i I 1 3 1 

(1) 

V fsdj 


where 


Tsd,. = 10^2 erg s"' 

(2) 

is the initial spin-down luminosity, and 


f,d = 2 X 10^ R-%B\IpI_^ s 

(3) 


is the spin-down timescale. Hereafter the convention (2x = 
2 /10^ is adopted in cgs units. 

The spin-down luminosity is essentially carried by a nearly 
isotropic Poynting-flux-dominated outflow. In the free zone 
(e.g. in the direction of the cavity drilled by the jet, or the 
intrinsically open regions in the ejecta), the magnetar wind 
would leak out from the ejecta, and undergo strong self¬ 
dissipation beyond R^ns > Rej > giving rise to extended emission 
(along the jet direc tion) or a bri ght X-ray afterglow emission 
(off-axis direction. IZhangll2013h . In the confined wind zone, 
the magnetar wind is expanding into the ejecta, and the mag¬ 
netic energy may be rapidly discharged via forced reconnec¬ 
tion (if R < /?dis) or self-dissipation (if R > Rdis)- The trapped 
dissipation photons would eventually show up when the ejecta 
becomes optically thin. 

As a rough estimation, one can assume an efficiency factor 
rji, to convert the spin-down luminosity to the observed lumi¬ 
nosity at frequency ly, so that one has 


4ttvD\ ’ 


(4) 


In this work, we take = 0.3 for both the extended emission 
and the late dissipated emission. 
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Fig. 1.— A cartoon picture of model, illustrating various emission sites and emission components at different stages, from a NS-NS merger event that results 
in a stable millisecond magnetar remnant. In the early stage, a short-lived jet was launched, punching through the ejecta shell and giving rise to the prompt short 
spike emission. Following the jet, the magnetar wind leaks out through the opening gap, dissipates at a larger radius, and powers the extended emission. In the 
intermediate stage, the magnetar spins down and the ram pressure overcomes the magnetic pressure, so that the gap is closed due to the hydrodynamical motion 
of the ejecta. The magnetar wind is trapped behind the ejecta, which heats and accelerates the ejecta, powering the merger-nova emission and the DNS merger 
afterglow emission. In the mean time, the initial jet energy would drive GRB afterglow emission. In the late stage, the ejecta becomes transparent, the magnetar 
wind still dissipates its energy and radiate X-ray photons, which freely escape the remnant and give rise to the rebrightening in the X-ray lightcurve. 


3.3. Magnetar powered merger-nova 

Suppose the magnetar is surrounded by a quasi-spherical 
ejecta shell with mass Mej and initial speed v,. A generic 
model for the dynamics and emission prope rties of the ejec ta 
could be briefly summarized as follows (e.g. lYii et al.l[2nill . 

Considering the energy injection from the magnetar and the 
energy dissipation through sweeping up the ambient medium, 
the total “effective kinetic energy” (total energy minus rest- 
mass energy) of the system can be expressed as 

£ = (r - 1 )Mejc2 + TE'^, + (r" - 1 )Mswc2 , (5) 

where F is the bulk Lorentz factor of the ejecta, is the in¬ 
ternal energy measured in the comoving rest frame, Mjw = 
^R^niUp is the swept mass from the interstellar medium 
(with density n) and R is the radius of the ejecta. The dy¬ 
namical evolution of the ejecta can be determined by 

dt MsjC^-\-El^^-\-2TMsv,c^ 

where T> = l/[r(l-/3)] is the Doppler factor with /3 = 
V The comoving time dt' and luminosity L' can be 

connected with the observer’s time and luminosity by dt' = 
Vdt and L' = respectively. 

With energy conservation, we have 

^ = ^L,, + VX^-V^L'^. (7) 


The radioactive power reads 


L;, = 4xl04‘'Mej,_2 


1 1 

- arctan 

2 TT 



1.3 

ergs-', (8) 


with fg ^ 1.3 s and ^ 0.11 s (iKorobkin et al.ll^l2h . The 
radiated bolometric luminosity L' reads " 




T > 1 


tR/T 
R/r ’ 


T < 1, 


(9) 


where r = k{M^^/V'){R/T) is th e optical depth of the 
ejecta with k be ing the opacity (iKasen & BildsterJ 120101 : 
iKotera et~ani2013h . 

The variation of the int ernal energy in the como ving frame 
can be expressed by (e.g. lKasen & Bildstenll2010l) 


rIF' HV' 

^ = CP-4,(10) 

where the radiation dominated pressure can be estimated as 
V = El^^/3V' and the comoving volume evolution can be fully 
addressed by 




( 11 ) 


The energy loss due to shock emission is ignored here, as is usually 
done in GRB afterglow modeling. 
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together with 


where the function g(p) takes the form 


dR 13c 
" (1-/3)' 


( 12 ) 


A full dynamical description of the system as well as the 
bolometric radiation luminosity can be easily obtained by 
solving above differential equations. Assuming a blackbody 
spectrum for the thermal emission of the mergernova, for a 
certain observational frequency ly, the observed flux can be 
calculated as 

1 2,tt^V^R^ {hvjVf 

^ 47 rD|max(r, 1) ]v’c'^v eKp(hv/'DkT')-V 


where h is the Planck constant. 


3.4. GRB afterglow emission 

The interaction between the initial launched jet and the am¬ 
bient medium could generate a strong external shock, where 
particles are believed to be accelerated, giving rise to broad¬ 
band synchrotron radiation (iGao et al.ll2013bl for a review). 
The total effective kinetic energy of the jet and the medium 
can be expressed as 

£ = (r-l)Mjetc2 + (r2-l)Mswc2, (14) 


where Ms* = 27r(l -cos 9)/3R^nmp with 6 being the half open¬ 
ing angle of the jet. The energy conservation law gives 


-(r^-i)(^) 

dt Mjet-l- 2 rMsw 


(15) 


where the energy loss due to shock emission is ignored. 

In the co-moving frame, synchrotron radiation power at fre¬ 
quency ly' from electrons is given by (Rybicki & Lightman 
1979) 


P', = 






(16) 


where is electron charge, = 3j^qeB'/{ATTm^c) is the 
characteristic frequency of an electron with Lorentz factor 7 ,., 
B' is the comoving magnetic field strength and 


F(x)=x / Ks/i(k)dk, 


(17) 


with K^i^(k) being the Bessel function. 

The comoving magnetic field strength B' could be estimated 
as 

B'= (18) 


where gj is the energy density in the shocked region and cb 
is the fraction of the shock energy density that goes into the 
magnetic field. 

The distribution of the shock-accelerated electrons behind 
the blast wave is usually assumed to be a power-law function 
of electron energy, 

dN' 

^0C7/, (7e.m < 7e < 7e.M), (19) 

Assuming that a constant fraction of the shock energy 
is distributed to electrons, the minimum injected electron 
Lorentz factor can be estimated as 

111 . 

7e.m = ^(p)e.(r-l)^, (20) 

me 


gip) 


P > 2; 


P-2 

In (7e.M/7e,m),P = 2. 


( 21 ) 


The maximum electron Lorentz factor 7 ^ m could be estimated 
by balancing the acceleration time scale and the dynamical 
time scale, i.e. 


7e.M 


TtqeB 
CnipC ’ 


( 22 ) 


where ^ 1 is a parameter that describes the details of ac¬ 
celeration. If the electron energy has a harder spectral index 
1 <p<2^ the minimum electron Lorentz factor w ould be de¬ 
rived as ([D ai & ChengllMTHBhattach'^l2()()!h 

_(2-pmp P- 2 V 

7e,m — ( . ee(r l)7e,M ) (^3) 

\p-l me ’ J 


With the dynamical description of the jet and these radiation 
equations, one can calculate the evolution of P',. Assuming 
that this power is radiated isotropically, then the observed flux 
density at frequency v = Vv' can be calculated as 


Fu = 


AttDI 


P', 

v'-i 


(24) 


Note that in the follo wing calculation, we neglect th e sideway 
expansion of the jet (lRhoadsl|1999HSari et ai.lll9'9^ . but con- 
sider the jet break with edge effect at later time when L"* > 0 
(iPanaitescu et aljri998ll . 


3.5. DNS merger afterglow emission 

During the propagation of the ejecta, a strong external 
shock would also form upon interaction with the ambient 
medium. With dynamical solution for the ejecta, and radi¬ 
ation equations [TS] to |24l one can easily calculat e the rele¬ 
vant br oad-band DNS merger afterglow emissions (iGao et alJ 
I2013ah . 


4. APPLICATION TO GRB 080503 

Considering the extremely deep limit for the host galaxy of 
GRB 080503, one plausible possibility is that GRB 080503 is 
at a moderately high redshift in a very underlumi nous galaxy, 
e.g. z « 1, comparable to the highest-z SGRBs (iPerlev et alJ 
l2009l) . In the following, we adopt z = 1 , and investigate the 
broad-band data of GRB 080503 with the physical model pro¬ 
posed in the last section. We And that with all standard pa¬ 
rameter values, the broadband data of GRB 080503 could be 
well explained. The fitting results (data from R. Hascoet et 
al. 2012) are presented in Figure |2] and the adopted parame¬ 
ter values are collected in Table 1. We briefly summarize the 
investigation results as follows: 

• The soft extended emission and the late X-ray excess 
peak could be well connected with a magnetar spin- 
down luminosity evolution function, suggesting direct 
magnetic dissipation as the same underlying origin for 
these two observed components. 

• In the X-ray band, the contribution from the merger- 
nova and early GRB afterglow emission are outshone 
by the aforementioned direct magnetic dissipation com¬ 
ponent; 
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Fig. 2.— Modeling results for the broad-band observations of GRB 
080503. The data are taken from R. Hascoet et al. (2012), with blue denoting 
X-rays and red denoting optical. Blue and red dashed lines represent the GRB 
afterglow emission in the X-ray and optical bands, respectively; blue and red 
dotted dash lines represent the merger-nova emission in the X-ray and optical 
band, respectively; the green dashed line denotes the evolution function of the 
magnetar spin-down radiation luminosity; the light blue dashed line denotes 
the late magnetar wind dissipation emission. The blue and red solid lines are 
final fitting lines for the X-ray and optical data, respectively. 


• The early optical data can be explained by the GRB af¬ 
terglow emission. The late optical data, including the 
re-brightening phase and the rapid decay feature, can 
be well explained by the emission from a magnetar- 
powered merger-nova. 

• Both the late-time optical and X-ray data peak around 
the same time when r = 1 , which is consistent with the 
argument that the late magnetar dissipation photons can 
travel freely after the ejecta becomes transparent. This 
powers the late X-ray excess; 

• For this particular event, the DNS merger afterglow 
emission is completely suppressed, since the ejecta 
mass is relatively large so that the ejecta is only mildly 
relativistic, and since the medium density is small. This 
emission component is not plotted in FiglJ] 

In the interpretation, we adopt the isotropic kinetic energy 
of the jet as 10 ^* erg, which is based on the total emission 
energy of the short spike and assume a factor of 20 % for 
the 7 -ray emission efficiency. The values for initial Lorentz 
factor (Fo) and half opening angle {9) of the jet are chosen 
as 200 and 0 . 1 , the values of which barely affect the final 
fitting results. To achieve the faintness of the GRB after¬ 
glow emission, a relatively low value for ambient medium 
density (n = 0.001 cm“^) is required, suggesting that GRB 
080503 may have a large offset relative to the center of its 
host galaxy, which in turn explains the extremely deep limit 
on its host at the afterglow location (since the system may 
have been kicked out far away from a host galaxy). The mi¬ 
crophysics shock parameters (e.g., e^, e^, C P) all cho¬ 
sen as their commonly used values in GRB afterglow mod¬ 
eling (iKumar & Zhangir20l4 for a review). For the magne¬ 
tar, a relatively large stellar radius /?s ,6 = 1.2 is adopted by 
considering a rapidly rotating supra-massive NS, and initial 


TABLE 1 

Parameters for interpreting the broadband data of 
GRB 080503, BY assuming z=1. 


Magnetar and ejecta parameters 

U(G) 

Pi (ms) Rs (cm) (M©) 

Vile 

K (cm’ g ‘) 

6 X I0‘^ 

2 1.2x10*^ 3x10“’ 

0.2 

10 

Jet and ambient medium parameters 


E (erg) To n (cm ’) 

9 (rad) 



I0’‘ 200 0.001 

0.1 


Other parameters 

Eg 

P C 



O.OI 

0.001 2.3 1 

0.3 

0.3 


spin period P, is taken as 2 ms by considering a mil d angular 
mome ntum loss via strong gravitational radiation (iFan et al.l 
l2013ah . The dipolar magnetic field of strength B is adopt 
as 6 X 10'^ G, which is consistent with the suggested values 
by fit t ing the SGRBs X-r ay plateau feature (iRowlinson et alJ 
1201 3t ILU & Zhang||2014l) . For the ejecta, we take the stan¬ 
dard values of mass (M^ ^ 1 O“^M 0 ) and initial velocity 
(v, = 0 . 2 c), and a relatively large value for the effective opac¬ 
ity k = 10 cm^ g“*. The latter was suggested by recent works 
by considering the bound-bound, bound-free, and free-free 
transi tions of ions (iKasen et al.ll2013l: iTanaka & Hotokezakal 
l2013h . If the opacity has a smaller value (e.g. because 
of the intense neutrino emis sion from the porto-magnetar, 
iMetzger & FernandezI (l2014h l. the same data can be inter¬ 
preted by increasing the mass of the ejecta. For example, 
an equally good fit can be reached with Mgj ~ 10 “^Mq for 
K = I cm^ g“'. Finally, we assume that 30% of the wind en¬ 
ergy is deposited into the eject a, which is a nominal value 
sugge sted from previous works (IZhang & Yanll20ll]; lYu et alJ 
l2013h . 

5. CONCLUSION AND DISCUSSION 

Double neutron star mergers could leave behind a millisec¬ 
ond magnetar rather than a black hole. In this scenario, 
the spin-down of the magnetar provides an additional energy 
source in the merger remnant, which generates much richer 
EM signatures from the merger remnant system than the black 
hole scenario. In this work, we give a comprehensive descrip¬ 
tion for all the possible EM signals under the magnetar rem¬ 
nant scenario, invoking several emission sites to account for 
several emission components, i.e. the initially launched jet 
to produce the short spike in prompt emission; an external 
shock site for this jet component to account for part of the 
observed optical afterglow emission component; an magnetar 
wind internal dissipation site that accounts for the early soft 
extended emission, the high latitude tail emission, as well as 
the late X-ray re-brightening emission when the ejecta be¬ 
comes transparent; an isotropic ejecta site that generates a 
magnetar-powered merger-nova emission; and finally the site 
where the ejecta interacts with the medium and powers the 
DNS merger afterglow emission. 

We presnt the detailed numerical methods to calculate these 
emission components and apply the model to investigate the 
broadband observations of GRB 080503. We find that the ma- 
gentar remnant scenario could well interpret the multi-band 
data of GRB 080503, including the extended emission and its 
re-brightening features in both X-ray and optical bands. In 
our calculation, we adopt z = 1 for GRB 080503, which could 
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be a plausible assumption in view of both the extremely deep 
upper limit for the host galaxy flux and the observed redshift 
distribution of SGRBs. If our interpretation is correct, some 
important implications could be inferred; 

• GRB 080503 is of a double neutron star merger origin; 

• The post-merger remnant of this event is a stable mag- 
netar, with an effectively polar cap dipole magnetic 
field 6 X 10'^ G and an initial period 2 ms; 

• The late optical re-brightening is a magnetar-powered 
merger-nova. Since its emission is essentially isotropic, 
similar merger-novae are expected to be associated with 
NS-NS merger gravitational wave sources even without 
a short GRB association; 


• For this event, the ejected mass during the merger is 
estimated to be around 3 x 10“^ Mq. 

To justify the assumption of z = 1, we also tested other red- 
shift values (either smaller or larger than 1). We find that 
the fitting results are not sensitive to the redshift value, even 
though some parameters may vary within reasonable ranges. 

In this work, we assume that the magetar wind is highly 
magnetized, i.e., with a high a value. If, on the other hand, the 
wind contains a significant fraction of primary pairs, the 
magnetic wind may become leptonic-matter-dominated upon 
interaction with the ejecta, so that a strong reverse shock can 
be de veloped, which would predict additional radiatio n signa¬ 
tures (lDail2004tlWang & Daill20f3tlWang et al.ll20l3l) . More¬ 
over, [MetzgeE^^^3( 20T3) proposed that the large optical 
depth of pairs inside the ejecta shell could also suppress 
the efficiency for depositing the wind energy into the ejecta, 
which essentially corresponds to a reduced value of ^ in our 
model. 

In our interpretation, we assume that the magnetar wind 
could leak out from the ejecta shell through the opening gap 
drilled by the initial jet, powering the extended emission. An 
alternative interpretation could be that the outflow from the 
magnetar wind itself may be collimated into a bipola r jet by its 
interaction with this ejecta (iBucciantini et al.ll2012l) and then 


power the extended emission (iMetzger & Pirol2014h . If this is 
the case, the real spin-down luminosity would be smaller than 
the extended emission luminosity due to the collimation ef¬ 
fect, inferring a somewhat lower dipole field. However, such 
a collimation effect is only significant for a large ejecta mass 
(say > 10“^ M 0 ), which should not affect the results in this 
work, since the preferred ejecta mass for the case of GRB 
080503 is relatively small (^ 10“^ Mq). 

Finally, it is worth pointing out that our described physical 
picture for the EM signatures from a NS-NS merger with a 
stable or supra-massive millisecond magnetar remnant could 
be applied to other cases of short GRBs and also the cases 
when the jet direction beams away from Earth. A system¬ 
atic study of extende d emission and internal plateau emis¬ 
sion from short GRBs (iLii et al.l2015h revealed many plateaus 
followed by a rapid decay. It would be interesting to sys¬ 
tematically apply the model to these GRBs to constrain the 
model parameters. In most cases, no X-ray rebrightening is 
observed, which suggests that the magnetar is likely supra- 
massive, and has collapsed into a black hole before the ejecta 
becomes transparent. In the future, off-axis X-ray transients 
may be discovered to be assoc iated with gra vitational wave 
events due to NS-NS mergers (IZhangll2013h . Applying our 
model to these events can give more detailed predictions to 
the brightness of these X-ray transients, which is valuable for 
searching for EM counterparts of GW signals in the Advanced 
LIGO/Virgo era. 
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APPENDIX 

The initial jet launched during the early accretion phase that powers the short GRB may have drilled a bipolar cavity in the 
ejecta. The subsequent magnetar wind following the short GRB also penetrate through this cavity power the extended emission. 
During this phase, the ram pressure around the cavity due to the dynamical motion of the ejecta would be initially balanced by 
the transverse magnetic pressure in the magnetar wind, i.e. 

^sd,i0 ... 

where c is the speed of light, (j) is the transverse magnetic pressure fraction, R is the radius of the ejecta, and A is the thickness 
of the ejecta shell. Under this condition, the corresponding fluid speed in the transverse direction due to dynamical motion of the 
ejecta can be estimated as 

/ 4d,0A y/" 

V A^ejC J 

^ 0.04c , (2) 

where A = v, Af, with v,- ~ 10'° cm s“' and Af = 1 s (IBucciantini et al.ll2012h . When t > t^d, the magnetic pressure quickly drops, 
so that the cavity would be gradually closed in a timescale 

^ Psd^sd^^ 

^close ~ 5 


V- 


(3) 
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where 6 is the jet opening angle, /3sd is the ejecta radial speed at fsd. From energy conservation, we obtain 


/3sd = min 


1 , 


^^rot 

MejC^ 


1/2 


For the cases with /3sd not close to unity (such as the case for GRB 080503 with ^ 0.5), we have 

feiose ~ 1.5 X 10'* s 


The total time scale for the extended emission can be estimated as 


(4) 

(5) 


tee ~ tsd + tclose. (6) 

With the parameters adopted to interpret the data of GRB 080503, we have fee ~ 74.4+ 157.3((^/0.25)“'/^ s, which is well 
consistent with the stopping time of extended emission (232 s), provided that the cj)^ 1/4 of the magnetic pressure concentrates 
in the transverse direction. 
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